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Abstract 

We calculate the thermally averaged rates for the r/-7r conversion and r] scat¬ 
tering using the Di Vecchia-Veneziano model and t’ Hooft model, which incor¬ 
porate explicitly the U(1)a anomaly. Assuming an exponential suppression of 
the U(1)a anomaly, we also take into account the partial restoration of U(1)a 
symmetry at high temperatures. We find that the chemical equilibrium be¬ 
tween rj and n breaks up considerably earlier than the thermal equilibrium. 
Two distinct scenarios for the 77 freeze-out are discussed and the correspond¬ 
ing chemical potentials are calculated. We predict an enhancement of the 
thermal ^-production as a possible signal of the partial U(1)a restoration in 

high-energy heavy-ion collisions. 
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I. INTRODUCTION 


At the Lagrangian level, QCD has, in addition to SU(Nf) x SU(Nf) chiral symmetry, 
an approximate U(1)a symmetry, under which all left-handed quark fields are rotated by a 
common phase while the right-handed quark fields are rotated by an opposite phase. It is well 
known that the I/(1)a symmetry is violated by the axial anomaly present at the quantum 
level and thus cannot give rise to the Goldstone boson which would occur when U{Nf) x 
U(Nf) chiral symmetry is spontaneously broken. The 17(1 )a particle, known as 7 /( 958 ) in the 
Nf = 3 case, acquires an additional mass through the quantum tunneling effects mediated 
by instantons [0], breaking up the mass degeneracy with pions and 7 / in the chiral limit when 
all quarks (u, d and s ) are massless. The 7/(547) particle also acquires an additional mass 
through the mixing with 7/. It is believed that at high temperatures the instanton effects 
are suppressed due to the Debye-type screening [0 . Then one expects a practical restoration 
of 17(1 )a at high temperatures. If the restoration occurs at a temperature lower than the 
chiral phase transition temperature T x , there may be some interesting phenomenological 
implications in high-energy heavy-ion collisions, as suggested first by Pisarski and Wilczek 
0 and more recently by Shuryak Q. One of the consequences of 17(1 )a restoration is the 
enhancement of 7/ particle production at small and intermediate transverse momenta due to 
the softening of its mass at high temperatures. However, the final yield of the 7/ particles 
and their p t distributions both depend crucially on the chemical and thermal equilibrating 
processes involving the 7/. 

In this paper, we shall examine the rates of various processes relevant for the thermal 7/ 
particle production, in particular, whether or not the 7/ can decouple early enough from the 
thermal system expected to be produced in relativistic heavy ion collisions. We shall present 
a theoretical calculation of the thermal cross sections for the processes 7/7/ 7/7/, 7T7/ th ] 

and 7/7/ 7T7T, essential to the thermal and chemical equilibration. Our calculations are based 

on models which explicitly incorporate the 17(1) a anomaly. We also assume an exponential 
suppression of the U(1)a anomaly due to the Debye-type screening of the instanton effect 
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0, which leads to the temperature dependence of the rj and rf masses. Our results suggest 
that the chemical equilibrium breaks up for rj particles long before the thermal freeze-out. 


We suggest a modest enhancement of thermal rj production as a signal for the relic of U(1)a 
restoration. 

This paper is organized as follows: In Sec. [ 11 ], we compute the mass spectrum of r/ and 
rf using the Di Vecchia-Veneziano model, which incorporates the U{1)a anomaly and the 
rj-rf mixing effect. We obtain the low-energy theorems for various scattering amplitudes. In 


Sec. Ill], we incorporate the a and the S resonances using the t’ Hooft model and reevaluate 


the r] scattering cross sections. In Sec. |TV], we study the thermal averaged cross sections 
responsible for maintaining thermal and chemical equilibria, and suggest that the chemical 
equilibrium between 7] and tt breaks up considerably earlier than the 7] thermal equilibrium. 
We discuss two scenarios for the 7] freeze-out and their corresponding signals for the tj 
production. We briefly comment on the roles of rf and the QCD sphalerons in Sec. 0 and 
Sec. [VI] respectively. 


II. NONLINEAR a-MODEL: LOW-ENERGY THEOREMS 

Up to now, there has been no direct experimental measurement of the 7] scattering cross 
sections (or the scattering lengths). One has to rely on theoretical models to calculate the 
interaction rates which are complicated by many uncertainties. Nevertheless, the scattering 
amplitudes at low energy can be more or less precisely predicted if the meson masses are 
soft, thanks to the soft-meson theorems which are based on the symmetry of the interactions 
and depend very little on the detailed dynamics. The current-algebra predictions of these 
scattering amplitudes have been made very early by Osborn || based on 577(3) x 51/(3), 
where the anomalous U{1)a and the r] and rj mixing are not included. In the light of 
softening of r] and rf masses at high temperatures, we argue that the symmetry can be 
extended to U( 3) x U( 3). We shall rederive the low-energy amplitudes incorporating the 
anomalous U{1)a using the nonlinear a-model that at the lowest order should give us the 
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low-energy theorems. The standard Di Vecchia-Veneziano model [6j.|7]], which incorporates 
the explicit U(1)a anomaly, reads after integrating out the gluon field 


n 2 p2 -P‘2 

£ cff = ^Tr (d^Udurf) + ^Tr (MU + MU ] ) + ^^-(iogdet U - logdetf/*) 2 , 
where U = exp (i$/f n ), f n = 93 MeV, M = diag (ml, ml, 2 m 2 K — m 2 ) and 


<f> 


1 7T° + Vv 7 3 + V2r]i/V3 V2n + y/2 K+ ^ 

\[2'k~ — 7T° + 7]s/V3 + V2 Vl /V3 V2K° 

v V2K- V2K° -2r/ 8 /V3 + V^rj^VS 


( 1 ) 


( 2 ) 


The last term in Eq. (© is the anomaly term which breaks U(1)a explicitly. It is easy to 
check that Eq. (]T|) satisfies the anomalous Ward identity which is crucial for determining 
the form of U(1)a breaking 0. In Eq. ([[]), a is related to the topological charge correlation 
function in pure Yang-Mills theory 


a = —i 


6 
J 7T 


d 4 x(T[F t , u F^(x)F liU F^( 0 )])ym , 


(3) 


where F^ v is the dual gluon field strength tensor and (• • •) stands for the vacuum expectation 
value at zero temperature or the thermal average at finite temperature. The integral a 
is identically zero in perturbation theory; it only receives nonperturbative contributions 
arising from the topologically nontrivial instanton configurations. The calculation of a at 
both zero and finite temperature is done by Gross, Pisarski and Yaffe [2j using a dilute gas 
approximation, and by Dyakonov and Petrov and by Shuryak [;9J using an instanton liquid 
model. For our purpose, the phenomenological value of a at T = 0 can be fixed by the meson 
mass spectroscopy, while a(T ^ 0) will be modeled by assuming an exponential suppression 
shown by Pisarski and Yaffe [ilQ] at high T. 

The quadratic terms for the octet rj § and the singlet rji from the Lagrangian reads 


r = 

A - , mass 


m 2 4 m 2 K ' 

IT + 3 


vl + 


'2m 2 K 


ml \ 9 2\/2 / 9 9 \ 

+ ^r + a )vi + ( 2 m n - 2 mr K ) mm 


■ (4) 


Clearly, there is a mixing between the octet rj 8 and the singlet r/i. The physical 77 (547) and 
77 '(958) are defined by 
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rj = 7]s cos 9 + rji sin 9 


r\ = —r] 8 sin 9 + 771 cos 9 


(5) 


to diagonalize the quadratic terms with the mixing angle 

. q 4 m K ~ ™l - 3 m 2 

tan 9 = -—— -- , 

2 \/2{m 2 K - ml) 

and the physical masses are 

m l — i m K + a / 2 ) _ \\f ( 2m |- _ 2 ml - a/3) 2 + 8a 2 /9 , 
m l r — ( m K + a/2) + i y/(2 m 2 K - 2ml - a/3) 2 + 8a 2 /9 . 


( 6 ) 

(7) 

( 8 ) 


The mixing angle 9 as well as ml and ml, depend on the instanton-induced quantity a which 


is a function of temperature. The precise form of a(T) at low temperature is not known. 
Nevertheless, if the U{1)a breaking becomes soft at a temperature lower than the chiral 
phase transition temperature T x , one may model the suppression effect by an exponential 
dependence [[T 0 | [T 2 | 


a(T) = a{ 0)e - (T/To)2 , 


(9) 


where T 0 ~ 100 — 200 MeV, while keeping the masses of the pion and kaon approximately 
temperature independent, since they change very slowly with the temperature. It is known 
that mixing angle 9, m^ and my at T = 0 cannot be simultaneously fit to their experimental 
values by a single parameter a(0). The best fit is to use the measured value of m^ + as 
an input to determine a( 0 ) = (mf ; + m^,) — 2 m 2 K and use this a( 0 ) to predict 9 , m^ and 
using Eqs. (§) and At T = 0, the predicted values are 9 = 18.3°, m r] = 500 MeV and 
rn v > = 984 MeV, compared to the measured values 9 exp ~ 20° from 77 , r / —> 77 , m exp = 547 
MeV and m e J p = 958 MeV. The temperature dependence of m v and m^ is completely 
determined by a temperature-dependent a(T) given in Eq. (^). Throughout this paper, 
we take To = 150 MeV in Eq. (0. It should be emphasized that in relativistic heavy-ion 
collisions, the thermal system freezes out at about T t h = 130 — 150 MeV, when the collision 
time scale exceeds the size of the system mainly determined by the nuclear radius R — 4 — 8 
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fm for central S + S or Pb + Pb collisions. Below the freeze-out temperature T t h, the hnite- 
temperature calculation of a(T) does not make sense and the behavior of a is determined by 
the nonequilibrium dynamics. Figure 1 schematically plots such a temperature dependence. 
Clearly, the 77 becomes soft at high T and eventually is degenerate with the pions. The mass 
of 7 / also decreases. However, it does not become degenerate with the pion because of the 
large strange-quark mass, as is seen from Eq. From Fig. 1 we see that the 7/ mass at 
high temperatures is still higher than the 77 mass at zero temperature. I 11 Fig. 1 we also plot 
the temperature dependence of the 5 resonance mass which we will discuss in the following 
section. At temperatures higher than T x , the masses of these excitation modes will all rise 
again. 

The interaction terms are obtained by expanding U in Eq. ([[]). In contrast to the pion 
held, there are no derivative couplings involving 77 and 7/. We shall ignore the interactions 
of 77 and 77 ' with kaons since they are heavy compared with pions. To the lowest order, the 
quartic terms involving tt, 77 and rf are 


A ini. — 


2 x 4 \fl 


2777,^.(77 sin y + v ' cos Xj + (4tt t 2 k — 2 ml) ( 77 ' sin y — 77 cos y) 


+ 12 m^ 7 r 2 (? 7 sin y + 77 ' cos 


( 10 ) 


where y = 0 + arctan l/\/2. At very high T, as 6 —* arctan \[2 and y —> tt/2, we can see that 
the 77 ' decouples from interactions with tt and 77 . The low-energy theorems on the two-body 
scattering amplitudes can be easily derived from Eq. (|T0|): 


^ 4.(7777 7777 ) 

4(7777 <->■ TT a TT a ) 

4.( 77 V 7T a 7T a ) 


f 

J 1 


2 [m v sin y + ( 4777 K - 2m n ) cos y] , 


1 


4(7T a 77 ^ TT a r]) = —7Ti n sin y , 
J 7 r 


4(7rV <-»• 7T a 7j') = —ml cos 2 y , 
J 7r 


4 ( 7777 ' 7T a 7T a ) 

rfrf) 

A(7]'r]' <->• 7777 ) 


4 (7r a 77 / ^ n a 7j) = — ml sin y cos y , 

J 7r 

^[tti 2 cos 4 y + (4m 2 K - 2m 2 n ) sin 4 y] , 

J IT 

4 ( 7777 ' 7777 ' = —(4 m 2 K — ml) sin 2 y cos 2 y , 

J 7 r 
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( 11 ) 


A{j]r] <r+ rjr}') 
A(rjrj' rfrf) 


— [ml sin 3 x cos x ~ (4 m 2 K - 2 ml) sin x cos 3 x] , 

H 

~p[ m l sin x cos 3 x ~ (4 m 2 K - 2 ml) sin 3 x cos x] • 

J 7T 


The results calculated by Osborn |>] based on the current algebra can be recovered by 
taking 6 = 0 and using the Gell-Mann-Okubo relation 4 m 2 K = 3m ; 2 ; + ml- These low-energy 
theorems must be satisfied by any dynamical models, because they are solely based on the 
symmetry properties of the theory. 


III. LINEAR (j-MODEL: INCLUSION OF RESONANCES 

The amplitudes listed in Eq. ( 0 ) grossly underestimate the strength of scatterings at 
higher energies, especially in the resonance regions. However, the inclusion of resonances 
introduces many uncertainties, such as which resonances should be included and what are 
the couplings of these resonances to the mesons. In addition, there is no guarantee that a 
naive lowest-order calculation will preserve the unitarity because of the strong interactions. 
Fortunately, the low-energy theorems provide us some guidelines as to how the amplitudes 
should approach their low-energy limits. The linear cr-model based on the chiral symmetry 
is known to satisfy the low-energy theorems, and at the same time to be able to incorporate 
the resonances. To further reduce the input parameters, we consider the a and 5(980) [now 
called do(980)] resonances, which, together with 7 r and the r] ns to be defined below, form a 
complete representation of U( 2 ) x 17(2). We shall concentrate on the 77 particle, since there 
is no dramatic change of the r( mass with temperature, as shown in Fig. 1. We study the 
most relevant processes for the 77 production: 7777 7T7T, tu] 7ir], and 7777 7777 . In this 

case, U(3)xU (3) reduces to 17(2) x U (2) except for the mixing effects which we have already 
calculated. 

Let us introduce the nonstrange mode r] ns = (uu+dd) /%/2 and take m s to be heavy. Then 
77 ns is approximately a mass eigenstate, 7] ns = 77 sin x + r f cos x, whose mass is determined 
from Eq. (f|) to be m\ s ~ 2a/3 + 777);. At zero temperature, m ns ~ 709 MeV. We then define 
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the (2,2) representation multiplet of 17(2) x U (2) as 


1 1 

$ = g ( a + ^ns) + ^ + ' i7V ) ' T ' 


The most general 17(2) x 17(2) invariant potential is 


( 12 ) 


V„ = + -(A, - A 2 )(Tr*tt) 2 + A 2 Tr(4> , 4>) 2 


(13) 


and the mass term is 


V m = 


™lf* 


Tr(<f> t + <f>) , 


(14) 


where Ai, A 2 are dimensionless constants. The 17(1) ^-breaking term, consistent with the 


Ward identity, is introduced by t’ Hooft [13] as 


K, = -(dett* + det*) , 


(15) 


and the coefficient in V a is chosen such that it gives the correct mass for r/ ns . The mass 
spectrum can be derived from Eqs. (0), (0) and (PD by making a shift a —> + a: 


2 \ c 2 1 2 2 \ r2 1 2 

™ a = X lJn + m n ; rn5 = A 2 f n + ? «ns • 


(16) 


The decay widths are 


r„ = 


32t r 


r<5 = {[K - ( m ns + m n ) 2 ][(m 2 s - (m ns - m n ) 2 }} l/2 


( m 6 ~ ? «ns) 2 


(17) 

(18) 


167r/^m! 

At zero temperature, T CT ~ 1 GeV (if ~ 700 MeV) and T^ ~ 200 MeV. In principle, 
we should also take into account the temperature dependence of /, r and m a below T x . 
Here, we assume the chiral phase transition is very rapid after which f n and rn„ have 
very slow temperature dependences. Furthermore, due to the large width of the a, the 
slow temperature dependence of m a will not change our results significantly. Under such an 
assumption, the linear cr-model predicts also some softening of the 6 resonance as T increases, 
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because S is the chiral partner of 7r and acquires some mass from the U(1)a anomaly. The 
temperature dependence of ms is plotted in Fig. 1. 

The interaction terms are 


Ant = + ^ns + <5 2 + 7T 2 )^ + ^fO 2 + Vis + ^ ^ 


2\2 


+A 2 fAo8 + 7?ns7r) • 5 + '-fiaS + j; ns 7T) 2 + 22(8 X 7r) 2 . 


(19) 


The coupling constants Ai and A 2 can be obtained from the mass relations of Eq. ([L6]). It 
is worth pointing out that the above model should not be used to estimate the pion-pion 
scattering amplitude, because it does not include the important vector resonances such as 
p and Ai. However, since pp and up scatterings cannot go through J = 1 channel, they 
do not directly affect the interaction rates for p. Similarly, we have also neglected the p-p 
interaction. 

To calculate the scattering amplitudes at the lowest order, we have to remove a pole 
singularity encountered when a resonance appears in the s-channel. A naive introduction 
of Breit-Wigner resonance width will spoil the delicate cancellation between the contact 
interaction and the pole exchange at low energy, leading to the violation of the low-energy 
theorems. We adopt a minimal prescription to save the low-energy limit developed by 
Chanowitz and Gaillard [[14], making the following replacement 


A if? 


Ai + 


s — mi + im a T c 


Ai(l - iT a /m 0 ) 


s — rnt 


s — mi + im a Y 0 


( 20 ) 


The scattering amplitudes are calculated as follows: 


A(pp «-► pp) = sin 4 X-A(PnsVns <-*• hnshns) 
= sin 4 %Ai(l - iT a /m a ) 
t — mi 


mz 


+ 


+ 


s — mi + im a T c 


u — mz 


t — rn 2 a + im a T a ' u — rri 2 + im a T 0 ’ 
A(pp <-> 7l a 7T a ) = sill 2 xAiPnsPns ^ 7T a 7T a ) 


= sin 2 yAi(l - iY^/mA) 


s — mz 


s - ml + 
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A(r]7T a <-> TjTT a 


+ sin 2 xA 2 (1 - iTs/ms) 
= sin 2 xA(rj ns n a <-> 77 ns 7 r“) 
= sin 2 xAi(l - iT a /m a ) 


t — mi 


+ 


u — m: 


t — m 2 + u — m 2 + imsT$ 

t — m 2 


sin yA 2 (1 — iTg/ms) 


t — ml + im a T 
s — m 


<7-*- <7 
2 


+ 


u — m: 


s — u — mj + im^T^ 


( 21 ) 


The cross sections for these processes are readily calculated by integrating out the scattering 
angle in u and t, most conveniently in the CM frame: 


a = 


f |P3c 


/-i 


\A\ 2 dcos9 , 


327rs |pi c 

where / = ( 1 ) 1/2 for (non-)identical particles in the final state. 


( 22 ) 


IV. THERMAL PRODUCTION OF THE ij PARTICLE 


We are interested in the production of 77 from a thermal source. To learn about the 
thermal history of the 77 , one needs to calculate the thermal averaged cross sections for 
various reaction channels. Since we are only concerned with the qualitative picture, we 
assume throughout the rest of this paper Boltzmann distribution functions for thermalized 
7 r’s and 77 ’s and ignore the quantum Bose-Einstein enhancement. The thermal averaged cross 
section for i + j —> k + / is 


(Vij<Tij{T)) 


1 (3,m i ,mj)A'i(v^/r) 

8 T mfmjK2(m i /T)K 2 {m j /T) 


( 23 ) 


where A(s, m*, mj ) = [s — (m, + mj) 2 } [s — (m ?; — mj) 2 } and ^/sq is the reaction threshold. The 
reactions 7777 —> 7777 and 7777 —> 7777 determine the collision time scale responsible for main¬ 
taining the thermal equilibrium while 7777 —» 7 T 7 t is responsible for the chemical equilibrium 
between 7 r’s and 77 ’s. We define the time scales r t h er and r c h em as 


r ther = (vcr(rjrj -> 777; /))n„ + (77(7(7777 -> vrvr))^ + (77(7(7777 -»• 7777)^ , 
r chem = {v<r{rm -> 7T77))7l r? , ( 24 ) 
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respectively, where n 7 r and n v are the number densities for n and rj, and the summation over 
different pion states is understood. We have performed a numerical integration in Eq. (|23| ) 
and plotted r t her and r c hem as functions of the temperature in Figure 2. In the calculation, 
we have explicitly taken into account the temperature dependence of m n (T), m ns (T ) 

and T g(T) as calculated in Sections II and III. We take a typical value R — 6 frn for the 
transverse freeze-out radius of the system. We define the thermal and chemical freeze-out 
temperatures T t h and T c h respectively as r t her(^th) = R and r c hem(Tdi) = R- One hnds from 
Fig. 2 

T th ~ 139 MeV and T ch ~ 168 MeV , (25) 


which are the temperatures at which the thermal and chemical equilibria start to break up, 
respectively. It is worth noting that T t h is comparable to the decoupling temperature of the 
thermal pions. 

The result that T c h is considerably higher than T t h offers an interesting possibility to 
detect the suppression of the U{1)a anomaly effect at high temperatures caused by the 
Debye-type screening. At sufficiently high temperatures T > T c h, the 7/ rescattering and 
the n-r\ conversion are frequent so that the system possesses both thermal and chemical 
equilibria. As the system expands and the temperature falls into the range T t h < T < T c h, 
the 7 r-T] conversion process becomes slow and effectively is turned off; the system can no 
longer maintain the chemical equilibrium. There is an approximate conservation of the total 
number of 7/’s since neither rjr) —> rjr) or nr] —> ni] can change the total //-number. The 
number density of ?/ at the chemical break-up temperature T = T c ^ is determined by the 
mass of r] at such a temperature m 71 (T c h): 


j [m"n (^ch)j A:h] 


1 

2^2 


m v (T ch ) 2 T ch K 2 


m^Tch) 

T c h 


(26) 


and the momentum distribution is just the Boltzmann distribution with zero chemical po¬ 
tential. However, this is not the final particle distribution, because the thermal collisions 
can still alter the momentum distribution. Nevertheless, the total number N v given by 
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N n = Ti? 2 r c n,[m,(T c h),T ch ] 


(27) 


is conserved at any time r < r c since rjr] tttt is turned off. Here m^Tch) — 360 MeV and 
r c is the proper time when the temperature of the system reaches T c h. 

As the system cools down to T t h, the mass of 77 should tend to m v (T t h) — 413 MeV, 
according to Fig. 1. If the rate for increasing m r? (T) is comparable to the thermal collision 
rate, the 77 particle adiabatically relaxes to m v (T t h). In this case, which we shall call Scenario 
A, one expects a standard thermal distribution for 77 at the freeze-out temperature T t h with 
a mass m v (T t h). The total number conservation requires 77 to develop a chemical potential 
H > 0 such that (neglecting the transverse expansion) 

T d e^ ,Tt ^n v [m v (T th ), T th ] = r c n v [m v (T ch ), T ch ] , (28) 

where r d is the freeze-out time when T = T t ^. The momentum distribution function in the 
local comoving frame is 

m \/ m ^( T th)+p 2 

/(p) = e^ /Tth e T th . (29) 


The chemical potential [i is a function of temperature, whose value at freeze-out can 
be determined from Eq. (|28|) once T d /T c is known. We assume that pions dominate the 
energy-momentum tensor (in fact we explicitly checked the contribution from // and found it 
negligible) so that T d /T c can be estimated by solving the ideal 1+1 dimensional hydrodynamic 
equation 




(30) 


where e is the energy density and P is the pressure, for massive pions. We find that T d / t c ~ 
1.53, given T ch /T th = 1.21. Substituting the ratio back in Eq. (|28|), one finds = e ,1 ^ Tth = 
1.58. We thus predict that if there is a partial U{1)a restoration at high temperatures, the 
thermal rj production given by Eq. ( |29| ) will be enhanced in this scenario due to both the 
finite chemical potential X v ~ 1.58 and a smaller 77 mass m 71 (T t h) — 413 MeV at the thermal 
freeze-out temperature T th . To quantify such an enhancement, we use Eq. fl29[) to calculate 
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the p t distribution of p particle, employing the fireball model and taking into account the 
transverse flow effects as described in Ref. ||T5|: 



where p = tanh 1 (/3 t ) and (3 t = (3 s {r/R) a (with f3 s = 0.5, a = 2) is the transverse flow 


velocity profile |L5|]. To reduce the possible normalization ambiguity, we also calculate the 


Pi-distribution for pions at the same freeze-out temperature TR, taking into account only 
the dominant resonance decays, p —> 27r, and plot the ratio 


V_ _ dN v /p t dp t 

7T° dN^/ptdpt 


(32) 


as a function of p t in Figure 3. It should be noted that the thermal ratio is only relevant 
when p t is small. At very large p t , hard processes become important and the fireball model 


is no longer applicable. For comparison, we also plot the same ratio for a normal case in 
which the p particles freeze out at the same temperature TR but with the zero-temperature 
mass m v = 540 MeV. 

Another situation, which we shall call Scenario B, is that the rate for increasing m v (T) 
when T t h < T < T c h is considerably smaller than the thermal collision rate. In this case, 
things get more complicated because the screening process is out of equilibrium. The p 
number conservation still holds, but the momentum distribution is quite different from that 
in Scenario A. Roughly one may imagine that even though the temperature drops to T t h 
after the chemical breakup, m v will still have the value m^(T c h), in close analogy to a 
“quenching” situation. The number density at the thermal freeze-out temperature is then 
n^[m, ? (T c h), Tth], and the chemical potential is determined by 


T d e ic/T ' h n v [m„(T ch ), T, h ] = T c n n [m,(T ch ),T ch ] , 


(33) 


yielding = e M//Tth ~ 1.24. The momentum distribution function is 


m \J m ^( T ch)+P 2 

/(p) = e^ /Tth e Air 


( 34 ) 
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which predicts larger 77 enhancement at low p t than at high p t . We also plot the ratio p/ir 0 
based on this scenario in Fig. 3 . 

What happens after the thermal freeze-out? It is clear that there must exist some 
mechanism for the 77 to relax from the ‘temporary’ entity whose mass is either m v (T t h) or 
rn v (T ch ) to its true identity at zero temperature with m v = 540 MeV. A possible picture 
might be that the 77 particles still feel a negative potential in the fireball. The height of the 
potential barrier is determined by the mass difference Am = m v — m v {T c h). The 7 ] particles 
with p t smaller than A?n will be trapped in the potential well until the rarefaction wave 
reaches the center of the interaction volume. Such a picture has been suggested by Shuryak 


[JJ and is similar to the mechanism of cold kaon production [ 18 |. At this stage, we do not 
attempt to address this nonequilibrium issue, but just to remark that our calculation here 
may have underestimated the enhancement effect at small p t £Am ~ 100 — 200 MeV. 

Both Scenarios A and B predict an enhancement of the thermal 7 ] production in the 
light of a partial U(1)a symmetry restoration. It would be very interesting to test the 
idea experimentally by measuring the ratio 77/7r°, especially its p t dependence. Although 
preliminary data from WA 80 [IB] on 77/77° ratio in both central and peripheral S + Au 
collisions at the CERN SPS energy, as indicated in Fig. 3 , have shown such a trend of 
enhancement, one certainly needs better statistics in order to make a definite conclusion. A 
related matter is the enhanced dilepton pair production via 77 Dalitz decay 77 —» If 

the 77 production is enhanced about 3 times, as we have predicted, the observed dilepton 


enhancement with the invariant mass below 500 MeV at the CERN SPS | 17 | may be partially 
accounted for. 


V. ROLE OF THE 77 ' 

There should be also some enhancement of the ratio 77'/7r°, since the mass of 77' also 
decreases as the temperature increases. Moreover, since the couplings of 77' to 77 and n in 
our model become small and eventually goes to zero when U(1)a is completely restored, 77' 
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might decouple from the system earlier than 77. The decay rf —> mrr] can also enhance the 
//-production. However, in our model, we postulate that the U{1)a restoration occurs at a 
temperature below the chiral phase transition temperature. Therefore, the kaon mass ttik is 
large and the rf does not become very soft. At T = T c h, the rf mass rry is about 750 MeV. 
Even without the effect of chiral symmetry breaking, the large strange-quark mass can give 
rise to a large mixing between 77 and rj according to Eq. (f|). This mixing gives rf a mass 
m 2 , = 2 m 2 K — m 2 even if U(1)a is completely restored. This mass is significantly larger than 
the 77 mass at any temperature. Therefore, in the context of our model, the rf effects are 
only moderate. The ratio r/'/ tt° should never exceed that of 77/77°. 


VI. ROLE OF QCD SPHALERONS 


So far we have confined ourselves to the possible suppression of the instanton effects at 
finite temperature that causes the softening of the masses arising from the topological charge 
transitions. At very high temperatures, it is known that such transitions can occur without 
going through the instanton configurations. In fact, they are dominated by sphaleron-like 
transitions whose electroweak counterparts have been extensively studied in the literature 
0. It is pointed out by McLerran, Mottola and Shaposhnikov 0 that the rate of a 
QCD sphaleron transition should be estimated in analogy to the electroweak theory for 
temperatures above the symmetry-restoration, which may not be quite suppressed. In the 
range of temperatures discussed in this paper, the rate of the QCD sphaleron transition may 
be unimportant. A rough estimate by Giudice and Shaposhnikov 


is 


r?™ = . 


8 K . 

~3~ 


(35) 


where k is the strength of the transition. The characteristic time scale of the sphaleron 
transition is 


50 


r sph = (192 na s T) ~ — . 


(36) 
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There is some evidence for k to be 0(1) from lattice calculations |2lJ. Unless k is really big, 
greater than 10, the sphalerons should be decoupled from the system in the hadronic phase, 
where the instanton effect is most dominant. 
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Figure Captions 


Fig. 1 The temperature dependence of rri 7] , m,y, ms- The parameter in the exponential 
suppression of the instanton effect is taken to be T 0 = 150 MeV. 

Fig. 2 The characteristic time scales of the thermal and chemical equilibration for the 7] 
particle. 

Fig. 3 The predicted ratio 77 / 7 r° as a function of the transverse momentum p t in three 
scenarios as discussed in the text. The preliminary data from WA80 |l(J are also 
indicated. 
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